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Abstract. We discuss a nonlinear plasma theory for self-modulation of pulsar 
radio pulses. A nonlinear Schrodinger equation is derived for strong electro- 
magnetic waves propagating in an electron-positron plasma. The nonlinearities 
arising from wave intensity induced particle mass variation may excite the 
modular. i <>nal instability of circularly and linearly polarized pulsar radiation. 
The resulting wave envelopes can take the form of periodic wave trains or 
solitons. These nonlinear stationary waveforms may account for the formation 
of pulsar microstructures. 
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1. Introduction 


Pulsar radio amissions exhibit ulti.ashurt intensity variations within 
individual pulses on time scales ranging approximately from 1 ps to 1 ms 
(Cordes , 1979). Since these microstructures (or micropulses) are the most 
fundamental features of pulsar pulses, they undoubtedly provide important 
clues concerning the physics of pulsar radio emission and propagation. 

Models for pulsar microstructure formation can be divided into two general 
classes: beaming models and temporal models. In beaming models, the rapid 

intensity fluctuation is considered to be a geometrical phenomenon caused by the 


angular sweeping of a nonuniform pulsar pencil beam across the observer's 1. ne 
of sight. For example, Benford (1977) views pulsar microstructures as the 
sweep of the observer's line of sight across a series of radiating plasma 
filaments directed along the pulsar magnetic field lines. Ferguson (1981) 
proposes that quasi-periodic microstructure emission is due to bunches of par- 
ticles located in many periodically spaced emission regions, perhaps occupying 
excited locations in a standing plasma v/ave of very long wavelength. In temporal 
models, the intensity fluctuation is treated as the result of an actual modu- 
lation of the pulsar radiation beam. For example, Cheng (1981) suggests that 
small oscillations in pulsar polar cap temperature can lead to strong modulation 
of the outflowing ions and electron-positron fluxes. Harding and Tademaru (1981) 
treat the temporal modulation of pulsar pulses as they propagate through magneto- 
spheric regions of relativistic velocity shear. 

In this paper we discuss a nonlinear plasma theory which may account for 
temporal modulation of coherent pulsar radio pulses. ' We demonstrate that the 
nonlinearities arising from relativistic particle mass variation may excite the 
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self-modulational instability of strong electromagnetic waves in an electron- 
positron plasma. The case of a circularly polarized wave is first studied in 
detail, then the treatment is extended to a linearly polarized wave. In general 
pulsar radio emissions contain both linear and circular polarization components 
(Cordes, 1979). The linear polarization originates from intiinsic emission 
processes such as the two-stream instability (Ruderman and Sutherland, 1975; 

Asseo et al., 1983), whereas the circular polarization may be converted from 

t 

the linear polarization via the Faraday rotation effect (Cocke and Pacholczyk, 

1980) or the nonuniform magnetic field effect (Hodge, 1982) as the wave propa- 
gates through the pulsar magnetosphere. 

Strong electromagnetic waves capable of driving plasma particles to rela- 
tivistic energies have received much attention in connection with the low-frequency 
pulsar magnetic dipole radiation (e.g.. Max, 1973; Kennel et al. , 1979; Chian, 
1982). Recently Bodo et al. (1931) used the dispersion relation of a strong 
linearly polarized wave to study the propagation of coherent pulsar radio 
emissions in the polar cap region. The high radio brightness temperature 
(10^ t-> 19 31 K, Cordes, 1979) inferred from pulsar observations implies that 
high-frequency pulsar radio waves can certainly drive the magnetospheric charged 
pai nicies to relativistic velocities. The intensity of pulsar radio emission 
in the source region can be measured by a dimensionless, Lorentz-invariant, 
parameter 


v 


eE 


m we 
c 


5.585 


D_ 

1 n J 6f 


( 1 ) 
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where S = flux density in Jy, f = pulsar emission frequency in MHz, Af = 
emission bandwidth in MHz, D = pulsar distance in Kpc, 6 = size of emission 
region in 10® cm, and n = index of refraction. Numerical examples of the wave 
strength parameter . are given in Table 1 for two pulsars, PSR 0950 + 08 
(Cordes and Hankins, 1979) and PSR 1133 + 16 (Bartel and Hankins, 1982), 
whose microstructure properties have been well analyzed. In Table 1 the index 
of refraction is taken to be unity, the observed peak flux densities are used, 
and the size of emission region 6 is determined by the product ct , where c = speed 
of light and t = pulsewidth. We expect that, 'with the exception of the pulsar 
distance, there may be several order of magnitude variations in most of the para- 
meters in Table 1. Hence it is reasonable to expect that the intensity of pulsar 
radio emission in the source region lies in the range 10 2 < v < 10 3 . For such 
wave intensities, plasma particles may acquire weakly relativistic to moderately 
relativistic velocities. This is in contrast to the case of low-frequency pulsar 
magnetic dipole radiation, for which the particles can reach ultrarelativistic 
velocities with v > 10 11 (Gunn and Ostriker, 1971). 

The modulational instability of electromagnetic waves has been extensively 
studied in connection with laser-plasma interaction for an electron plasma (e.g.. 
Max et al., 1974; Shukla et al . , 1977) and an electron-ion plasma (e.g., Berezhiani 
et al., 1980). While these studies are surely of interest to laser fusion appli- 
cations, they might not be appropriate for pulsar applications. According to 
current polar-cap pulsar models (Ruderman and Sutherland, 1975; Arons and 
Scharlemann, 1979) the pulsar magnetosphere is composed of secondary electrons 
and positrons resulting from pair production induced by high energy curvature 
radiation photons emitted by primary positron or electron beams coming from 
the pulsar surface. In such a magnetosphere, positrons and electrons contribute 
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equally toward the wave characteristics (Clemmow, 1974; Kennel and Pellat, 1976; 
Chian, 1980) , hence the inclusion of positron dynamics is essential for pulsar 
radio emission theories. Moreover, the ion acoustic mode is absent in an electron 
positron plasma (Tsytovich and Wharton, 1978), therefore the oodulational coupling 
between the high-frequency electromagnetic wave and the low-frequency ion acoustic 
wave, heretofore considered for an electron-ion plasma, may not be applicable to 
the pulsar environment. 

Astrophysical applications of the modulational instability of electromagnetic 
waves in a hot plasma have been studied for the situation where the relativistic 
effect is caused by high plasma temperatures (Pataraya and Melikidze, 1980; 

Lakhma and Buti, 1981). In the present paper we adopt the model of a cold, 
unmagnetized plasma in order to single out the effect of relativistic particle 
mass variation originating from the interaction of electrons and positrons 
with high intensity waves. Although the magnetic field is certainly not neg- 
ligible in the pulsar magnetosphere, we postpone the analysis of its effects to 
another paper in order to investigate the simplest case in detail. 


2. Theory 

Consider the propagation of r. circularly polarized electromagnetic wave in 
the rest frame of an electron-positron plasma with its vector potential given by 


A(z,t) - a(z,t) (x cos tot + y sin tot) 


C2) 


In the linear regime the wave is governed by the dispersion relation 


to 2 = c 2 k 2 + to 2 

P 


tO 2 = 10 2 + l0 2 

p pe pp * 


( 3 ) 


where = *2 is the electron-positron plasma frequency, = w ’ = 

(4nN e 2 /m 
O 6 

The nonlinear wave propagation is described by the following wave equation 


— — - c 2 V 2 A = 4 Tree (N v - N V ) 
at 2 ~ P'P e ~ e 


(4) 


where only purely transverse fields are treated, since we may neglect the longi- 
tudinal electric field in an electron-positron plasma (Kennel and Pellat, 1976). 
The relativistic equations of motion for electrons and positrons are 


a q v ■< B 

(_ + y • V) (y v ) = — (E + — ) 

v 3t *a • v a~a m c 

e 


'5) 


where a = (e,p) and y = (1 - v q 2 /c 2 ) From (2) and (5) we have 


e A 

~ e ro e c(l + XA 2 )^ 


e A 


v = 


' p m.^cfl + XA 2 )' 


( 6 ) 


where X = e 2 /m e 2 c l *. Making use of the quasi-neutrality condition, N e ~ N p = N c » 
appropriate for a circularly polarized wave the nonlinear wave equation becomes 


3 2 A 3 2 A 

,2 1 


3t 


- c" 


+ to 


3z " P (1 + XA 2 )^ 


= 0 


(7) 


Introducing a complex modulational representation 


A(z,t) =j[a(z,t)e 10Jt + c.c.] 


( 8 ) 
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(7) can be rewritten in terms of the slowly time-varying modulation 
a (with a 1 3a/3t << w) 


da 2 3 a w z 

c a) p 

1 — + — TTT * — a - 


at 2w az z 


2 - 


2io 


(1 + X\q\ 2 f 2 


= 0 


( 9 ) 


where a term proportional to 3 2 a/9t 2 is dropped. It can be shown that (9) 
admits localized stationary solutions for arbitrary wave amplitudes, however, 

J 

for the sake of illustration we shall henceforth treat the small amplitude 
limit of (9) . 

In the limit l|a| 2 « 1, (9) yields a nonlinear Schrodinger equation 


3a 

i — + P 
9t 


3 2 a 

8z2 
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Q = 
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An additional term, Ra (where R - (u 2 - w p 2 )/2w), was removed from (10) by the 
transformation a a exp(iRt). Equation (10) is the desired wave equation for 
the modulational instability. Alternatively, (10) can be derived by 
a method introduced by Karpman and Krushkal (1969), by which the coefficients 
of the nonlinear Schrbdinger equation (10) are determined from a weakly 
nonlinear dispersion relation 


on = w(k, |a| 2 ) 


with 


P = 


1_ 3 2 o) 
2 3k2 


Q = - 


3 03 


3liP r 


( 12 ) 


(13) 




The exact nonlinear dispersion relation for a circularly polarized electro- 

_ ^ J _ />! n /«4- *>Ari niDCTna -»< if* o J.r»0»J • 
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2 V s (1 + v 2 )'* 


( 14 ) 


where v 2 = e 2 E 2 /m 2 io 2 c 2 = XA 2 . Note from (14) that the positron contribution 
e 

to the wave dispersion is equal to the electron contribution. In the weakly 
nonlinear limit v 2 << 1, (14) reduces to . 

w(k, |a | 2 ) = (c 2 k 2 + o> p 2 - \a\ 2 ) h C«) 

Applying (15) to (13) we obtain precisely the same coefficients as those in 
(11). It is important to observe that the Kaipman-Krushkal method is only 
valid for the weakly nonlinear regime; it should not be applied directly to 
a strongly nonlinear dispersion relation such as (14) (see e.g., Durrani et al., 
1980). The appropriate modulational wave equation for large wave amplitudes 

is (9). 

Having derived the nonlinear Schrodinger equation, (10), we can now determine 
whether or not the plane carrier wave (2) is unstable to a low-frequency modu- 
lation. It has been established (e.g., Hasegawa,, 1975) that the modulational 
instability can be excited in a medium provided the group dispersion P and 
the nonlinear frequency shift Q are of the same sign, namely PQ > 0. From (11) 

we see that this condition is indeed satisfied. 

Consider next the dynamics of the modulated wave moving with a group speed 
V = 3w/3k, described by (10). In this case we may assume a finite value a Q 
-or the modulation as z-Vt-*- ± thus in (10) \q\ 2 is replaced by |g| 2 - l<2 0 l 2 
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1 aT + P + Q( l~l 2 " = 0 * (16) 

The stability of the modulation can be studied by separating a into two real 
variables, p and a, representing the real and the imaginary parts of a 

a = /p (z,t) - (17) 


Substitution of this expression into (16) gives for the real and imaginary 
parts, respectively. 


QCp - & “f 


JL (IB ., 2 

4p2 W 
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3t 


= 0 


(18) 




(19) 


Linearizing (18) and (19) 


i (k. z - Ur t) 
p = p + p e v L L 

O 1 


( 20 ) 


a = 0l e i(k L z ' 


( 21 ) 


we obtain the following dispersion relation for the low-frequency modulation 


to. 2 = P 2 k. 4 - 2PQ-p k 2 
L L o L 


( 22 ) 
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It follows from (22) that the threshold for the modulational instability is 
P Q > Pk L 2 /2Q, namely 


v 

o 


> ck /w 
L p 


( 23 ) 


for the circular polarization case and the corresponding maximum growth rate 
is y = Qp q , namely 

Y = Up 2 v o 2 /4w . • (24) 

The instability is purely growing, and therefore non-convective in the rest frame. 

As the modulation grows, the instability can evolve to a nonlinear station- 
ary state that results from the balance between nonlinearity and dispersion. 

The possible stationary solutions are the periodic wave trains (Nishikawa a’.d 
Liu, 1976) 

P . p o cn2[|^-| % P 0 C2 - Vt)] , (25) 

or the solitary waves 

P = P Q sech 2 [|-^r|^P o (z - V.t)] . t 26 ) 

The solution, (26), is called an envelope soliton since it is the envelope of 
the wave packet that has the form of a solitary wave. These envelope solitons 
can be shown to be stable against longitudinal perturbations and mutual collisions 
We now apply the Karpman-Krushkal method to discuss the self-modulational 
Instability of nonlinear linearly polarized electromagnetic waves in an electron- 
positron plasma. Although the dispersion relation in the strongly nonlinear 
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limit has been derived by Kennel and Pellat (1976) and Bodo et al (1981) 
it cannot be used here, since the Karpman-Krushkal method is restricted to 
the weakly nonlinear regime. The nonlinear dispersion relation for a linearly 
polarized wave in an electron plasma, in the limit v 2 << 1, is given by 
(Chian and Cletmnow, 1975) 

_ to 2 

n 2 = 1 - (1 - v 2 ) -E£- . (27) 

8 o ) 2 

Now, (9) suggests that the dispersion relation for an electron-positron plasma 

can be obtained from the dispersion relation for an electron plasma by replacing 

id for cd (i.e., the electron-positron plasma frequency cd ) . Hence the 
pe pe P 

appropriate dispersion relation can be obtained from (27) by replacing (b for 

i 

cDp, namely 

3cd 2 X , 

otk.lal 2 ) - (c 2 k 2 * - — §— lal 2 )' 4 . - (28) 


A comparison of (28) with (15) sh.-ws that, apart from a slight difference in 
the numerical coefficient, the dispersion relations for linearly and circularly 
polarized waves are essentially the same. Applying (28) to (13) then gives 
a nonlinear Schrodinger equation (10) for a linearly polarized wave with 



Q = 


3 e 2 (D- 2 

E 

16tDm 2 c 2 
e 


(29) 


Evidently, PQ > 0, thus the self-modulational instability of linearly polarized 
waves in an electron-positron plasma can be excited. The equations describing 
the evolution of the modulational instability for this case are similar to the 
circular polarization case with P and Q given by (29). 



3. Discussion 

In the previous section we demonstrated that the self-modulational 
instability of circularly and linearly polarized electromagnetic waves can be 
excited in an electron-positron plasma. The theory presented in this paper is 
compatible with the statistical model for the pulsar signal proposed by Rickett 
(1975). In Rickett' s model the pulsar signal is depicted as an amplitude- 
modulation a (t) of a noise-like random process n(t): 

p(t) = a (t) n(t) , 

where p(t) is the electric field of the received pulsar signal, n(t) describes 
the coherent fast time-varying emission by particle bunches, and a(t) describes 
micropulses that vary much more slowly than n(t). This amplitude-modulated 
noise model has been found to be consistent with the observed pulsar spectra. 
According to our theory, as evidenced by (8), the high-frequency coherent 
pulsar emission with a fast-time scale 2tt/w can be modulated by a slowly time- 
varying envelope g(z,t) due to nonlinearities arising from relativistic particle 
mass variation. Hence, the resulting modulation envelopes may account for the 
formation of pulsar microst s/uctures . 

We have seen that the end products of the relativistic self-modulational 
instability may be nonlinear stationary waves of either periodic wave-train 
type or envelope soliton type." This is consistent with the observed features 
of pulsar microstructures. Hankins and Boriakoff (1978) showed that the observed 
pulsar microstructures can be put into two categories, intermittent type or 
quasi-periodic type, according to their intensity structures. Most micro 
pulses are of the intermittent type, which have bursts of strong emission 


interpersed with sections of longitude where the signal returns abruptly to 
the system noise level. Occasionally, quasi-peri odic string of micropulses 
are detected. The intermittent micropulses can be explained by our theory as 
a collection of envelope solitons with randomly fluctuating amplitudes, whereas 
the quasi-periodic micropulses can be explained as periodic wave-trains or a 
sequence of envelope solitons with little variation in their peak amplitudes. 

We now make use of the nonlinear envelope solutions 


A(z,t) 


A 0 cn 2 [|^-|^ A q (z - Vt)] 

< 

A o sech 2 C|^| % A o ( z - Vt)] 


(31) 


to estimate the features of microstructures produced by the relativistic self- 
modolational instability. First we calculate the Lorentz-invariant number N 
denoting the number of wave crests in a given modulation. From the observed 
temporal pulsewidths of pulsar micropulses 1 us to 1 ms) and the frequency 
range of put sar radio emission (~ 100 MHz to 1 GHz) we see th-'.t N varies roughly 
from 10 2 to 10 6 . (31) shows that the envelope spatial pulsevidth for the circular 
polarization case is 


i2P|Js _1_ _ 2c . 

■O' A ai v 
x o p o 


Thus, in the rest frame. 


N = - 2 — — = -i— (— ) 
0 ) v 2 tt 7TV V ( 1 ) j 
p o o p 


(32) 


( 33 ) 


where k is the wavenumber of high-frequency pulsar radio waves. (33) shows 

-i r> rnr W fn n p w i Hiifi flu* oi'i^f-TV^d TrtTiPG I l IS rCGuirCu that 

V. A** **-“ “ ~ ' W 

k'c/io >> 1. The weakly nonlinear dispersion relations (15) and 

^ . -2 
(28) indicate that this is possible only if w >> For example, if v = 10 

and a) /u = 10, N = 10 3 . Therefore, our theory suggests that self-modulational 
P 

formation of microstructures takes place in regions of pulsar magnetosphere 


where the pulsar radio wave frequency is considerably greater than the local 

plasma frequency. Next we use our theory to calculate the temporal pulsewidth 

of micropulses. Since (32) is given in the rest frame of secondary particles 

it is necessary to transform (32) to the observer's frame. Because of time 

dilation and Doppl ;r shift (Ruderman and Sutherland, 1975) in the observer's 

frame uj + 2y u . where y is the relativistic factor associated with the 
’ pe s pe s 

stream velocity vf secondary electrons and positrons. Hence, in the observer's 
frame, the temporal pulsewidth is 


t - 


y to v 
s pe o 


(34) 


where the relation 6 = ct is used, since for y » 1, the soliton moves with 

b 

relativistic speeds in the observer's frame. If = 10 2 , w pe = 10 6 s -1 (in the 
rest frame) and v = 10“ 2 , then t = 1 ps which is within the observed range. 

The above calculation suggests again that relativistic self-modulational insta- 


bility is excited in the region of magnetosphere with relatively low plasma 
densities. Now (34) shows that micropulses of higher intensity have narrower 
pulsewidths. Such behavior was shown by Bartel et al. (1980) to hold true 
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for pulsar subpulses and was suggested by Ferguson (1981) to be true also for 
micropulses. It would be interesting for pulsar observers to check this property 
across the entire range of microstructure time scales. 

In conclusion, we point out that the physical mechanism for formation of 
pulsar microstructures discussed in this paper is expected to take place in 
most pulsar magnetospheres since the wave intensity induced relativistic effect 
responsible for driving the self-modulational instability is general. However, 
a quantitative test of a modulational instability theory requires the extension 
of the simple model adopted in this paper to include other effects, such as 
large wave amplitude, ambient magnetic field, plasma temperature and plasma 
inhomogeneity . In the meantime, cur simplified analysis indicates that the modu- 
lational instability is a promising direction for further investigation. 
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TABLE 1 


Numerical examples of the wave strength parameter 


— — '"1 

Pulsar 

S(Jy) 

/(MHz) 

A/ (MHz) 

D(Kpc) 

6(10® cm) 

mm 

PSR0950 + 08 

850 

430 

200 

0.1 

6 x 10' 3 

89 

PSR1133 + 16 

300 

1720 

10 

0.2 

7.5 x IQ -4 

47 











